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We study the renormalization of the gap of an s-wave superconductor in the presence of two 
magnetic impurities. For weakly bound Shiba states, we analytically calculate the part of the 
gap renormalization that is sensitive to the relative orientation of the two impurity spins. For 
impurities with a strong exchange coupling to the conduction electrons, we solve the gap equation 
self-consistently by numerics and find that the sub-gap Shiba state turns into a supra-gap Andreev 
state when the local gap parameter changes sign under the impurities. 
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The study of magnetic impurities in superconductors 
has a long history. Important effects include the renor¬ 
malization (reduction) of the gap in the host supercon¬ 
ductor by the impurities^i^ and the induced Yu-Shiba- 
Rusinov, or Shiba, bound statesi^^— In recent years, re¬ 
newed interest in chains of magnetic impurities in super¬ 
conductors was furthermore driven by their potential of 
hosting Majorana zero modes.— 

Using scanning tunneling microscopy, the Shiba states 
induced by individual magnetic impurities on the surface 
of a superconductor have been shown to be strongly lo¬ 
calized at the impurity sitesi^d^ This can be understood 
as a consequence of their sub-gap energy. When this 
energy is below the Fermi level, the Shiba state is occu¬ 
pied by a single electron, whose spin direction is dictated 
by the impurity spiui^ Shiba states are to be contrasted 
with Andreev bound states in extended superconductor- 
normal-superconductor (S-N-S) heterostructures, or in 
heterostructures composed of superconductors with dif¬ 
ferent gaps (S-S’-S). Andreev bound states can be viewed 
as standing waves of electrons and holes inside the gap¬ 
less region that result from Andreev reflections on the 
S-N-S and S-S’-S interfaces.—^— 

In this work, we address a scenario intermediate be¬ 
tween single impurities and impurity chains, and specif¬ 
ically study the renormalization of the superconducting 
order parameter in the presence of two magnetic impu¬ 
rities. We quantify the inter-impurity scattering analyt¬ 
ically and discuss how the Shiba states are affected by 
the gap renormalization. In particular, we observe that 
sub-gap Shiba states can transmute into supra-gap An¬ 
dreev states when the impurities are strongly bound to 
the superconductor, which we show by a simple analyti¬ 
cal model and by self-consistent numerics. 

The outline of the paper is as follows. In the next sec¬ 
tion, we introduce our model Hamiltonian for the mag¬ 
netic impurity in an s-wave superconductor. In Sec. |lll 
we summarize briefly our analytical approach based on 
the T-matrix formalism. In Sec. imi we apply it to a sin¬ 
gle classical impurity. In Sec. IIVI we present a detailed 


study of the Shiba to Andreev transition as a function of 
the impurity strength using both analytical and numer¬ 
ical methods. In Sec. El we consider the two-impurity 
case, and show the presence of two subsequent quantum 
phase transitions as a function of the exchange coupling. 
Finally, in Sec. ED we present a summary of our results 
and some perspectives. 

I. MODEL HAMILTONIAN 

In the analytical part of our study, we analyze 
two (weakly bound) magnetic impurities in a three- 
dimensional s-wave superconductor. The corresponding 
Hamiltonian H = iJkin + Hqc + ^^imp can be decomposed 
into the kinetic energy, i/kin, the BCS mean field super¬ 
conductivity, i?sCi and the impurity contribution, Hmip- 
If there were no impurities, the system could be modeled 
by the Hamiltonian Hq = Hkin + Hgc.o, where Flkin = 
Ek <7 + ^.c. (Uk denotes the kinetic energy 

dispersion), and Hsc,o = / dr Aq cj(r)c|(r) -|- H.c. with 
a superconducting order parameter of constant value Aq 
( chosen to be positive). Here, Cka is the annihilation op¬ 
erator for an electron of three-dimensional momentum k 
and spin cr, while cj.(r) creates a spin a electron at posi¬ 
tion r (in the calculations, we use t= +, 1= —)■ Due to 
the impurities, however, the superconducting part of the 
Hamiltonian is promoted from idsc.o to 

Hsc = j dr A(r) c|(r)c|(r)-k H.c. (1) 

The order parameter is set by the self-consistent equation 
A(r) = —g {ci{r)c^{r))H, where 5 > 0 is the microscopic 
attractive interaction between electrons in the supercon¬ 
ductor, and where the expectation value is taken with 
respect to the full H. The magnetic impurities, finally, 
are treated as classical spins Si residing at positions ri, 
and modeled by a purely magnetio^ and point-like scat¬ 
tering potential. Their classical treatment implies that 
the we restrict our analysis to sufficiently large impurity 
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spins. We address impurities polarized along the z axis 
in spin space, thus covering both parallel and antiparallel 
impurity spin alignments. The corresponding impurity 
Hamiltonian reads 

^ ^ Ji^iz ^ (^2)^(T(^0 2 (^) 

2,(T 


space (tq = I 2 X 2 ), while we obtain 


So{r,r',uJn) 


nvF 


( 6 ) 


X ^sin(fcF<5r) (ia;„ro + AqTx) + cos{kFSr)^J^ 



where Siz = ±|Si| is the z component of the spin of 
impurity i, and Ji is the exchange coupling between this 
impurity and the electrons in the superconductor. This 
model comes with the Debye frequency wd as a natural 
high-energy cutoff. 

If the renormalization of the superconducting gap is 
small, |dA(r)| = |A(r) — Ao| Aq, one can approxi¬ 
mate A(r) by evaluating the right-hand side of the gap 
equation for an unrenormalized gap, 

A(r) Ri-g(c^(r)ct(r ))//2 (3) 

with H' = Hq + Hiznp- This is the approach we em¬ 
ploy in the remainder of the analytical calculation. When 
|dA(r)| becomes of order Aq, this approximation ceases 
to be valid, and we make use of numerical simulations. 
The tight-binding Hamiltonian H is defined as 

H = -t (4) 

<z,i'> cr —±1 

N^-Ny 

+ ^ [AjaCiaCia Jia)cl^Cia] + H.C., 

i=l a=±l 

where Cia is the annihilation operator acting on an elec¬ 
tron with spin a at lattice site i, and the first sum 
runs over neighboring sites i and i' located in a two- 
dimensional square lattice of size Nx x Ny with lattice 
constant a. The chemical potential ^ is taken from the 
bottom of the energy band, and the local order parameter 
Ai is determined self-consistently in an iterative proce¬ 
dure for fixed value of the exchange coupling Ji at the 
site i starting from the uniform superconducting order 
parameter Aq until convergence is reached^^— 

II. ANALYTICAL T-MATRIX APPROACH 


at distances 5r = |r—r'| larger than the Fermi wavelength 
(kF is the Fermi momentum). Using the imaginary time- 
dependent Dyson equation, the full Green’s functions can 
be expressed as 

g{v,r',uJn) = go{r,r',uJn) 

+ Goirj,r',u}„) , (7) 

ij 

where the T-matrix Tn = V [l 4 x 4 ~ Gp^nV] ^ involves 
the 4x4 -matrix Go,n with entries Gq^ = Wn), 

and V = dia,g(JiSiz, JiSiz, J 2 S 2 Z, J 2 S 2 z)- 

III. SINGLE IMPURITY PHYSICS 

When the impurities are far apart from each other, 
{kFri 2 )~^ ^ 1, the gap near a given impurity is pre¬ 
dominantly renormalized by scattering processes off this 
impurityiii^ We can then for instance focus on impurity 
1. Next to this impurity, and using ai = nvFJiSiz, the 
dominant contribution to the gap renormalization at ri 
reads^ 

( 8 ) 

ggfA q f°° + al) + JA^ 

2 i-00 ^ ■ 

A Gaussian high-energy cutoff at the scale ujjj regularizes 
the logarithmic UV divergence of this integral. For small 
Ioi I, we can furthermore expand this expression to second 
order in ai. Since the self-consistent equation ([3]) yields 
a bare gap of Aq « 2a;D we find 

dA(ri) R! -3a^Ao . (9) 


We begin with calculating the full imaginary 
time Nambu Green’s function C/(r, r', r, r') = 

— using the ‘bare’ Green’s 

function go(r,r',T,T') = -(T^^'(r',(r,r))//o, and 
where the Nambu spinor is defined as = (Ck^i'i 
For equal positions r = r', a Fourier transformation 
from imaginary time to Matsubara frequencies yields 

So(r,r,uJn) = - (iujnTo + AqTx) , (5) 

V^n + 

with vf being the density of states (per spin) at the Fermi 
energy, and with the Pauli matrices Ti acting in Nambu 


For a more strongly bound Shiba state, |q;i| —?► 1, our an¬ 
alytical calculation hints at a very strong suppression of 
the gap at the positions of the impurities, |( 5 A(ri)| ^ Aq, 
but the non-selfconsistent approach cannot make any 
quantitative predictions in this regime. Numerically, on 
the other hand, we can readily access the regime |q;i| > 1 , 
see Fig.[T] For small values of Ji, the gap under a (single) 
impurity is suppressed quadratically in Ji in agreement 
with Eq. (|n]). Interestingly, at a critical value Jc corre¬ 
sponding to the phase transition,—!^ the superconducting 
gap under the impurity changes its sign and magnitude 
abruptly, giving rise to a local tt junction, see Fig. [I}i. 
Also the energy of the bound state jumps from zero to 
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a finite negative value. This feature results from self- 
consistency in determining A(r) and cannot be captured 
analytically in the approach we used above. The stronger 
Ji becomes, the smaller Ai under the impurity gets, until 
the superconducting gap is eventualy totally suppressed. 
As a result, the internal S-S’ junction evolves into an S-N 
junction. 
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IV. SHIBA TO ANDREEV STATE TRANSITION 


After the phase transition, the modulus of the bound 
state energy increases, whereas the superconducting gap 
I All decreases as functions of Ji. This means that the en¬ 
ergy of the bound state lies outside the gap for Ji > Jsa- 
We refer to this state as an Andreev bound state local¬ 
ized in an effective S-S’ junction. The transition is clearly 
visible in our numerical results plotted in Fig. [1^. A key 
signature distinguishing Shiba and Andreev bound states 
is their degree of localization. This is demonstrated in 
Fig.l which shows the numerically obtained bound state 
wave functions. Indeed, the Shiba states are more local¬ 
ized under the impurity, see Fig. [2^. In contrast, the 
Andreev bound states are delocalized directly under the 
impurity as their energies are higher than the supercon¬ 
ducting gap in this region, see Fig. This difference 
in the characteristics of the bound state wavefunctions 
can be checked experimentally with STM techniques. We 
note that the hlling of the Shiba state does not necessarily 
coincide with the Shiba to Andreev transition, see again 
Fig. [T^. As a consequence, there exists a finite range of 
exchange couplings with a filled but localized Shiba state, 
which delocalizes at the Shiba to Andreev transition. Our 
results agree qualitatively with earlier numerical studies 
of a single impurity^^d^ which, however, do not address 
the Shiba to Andreev state transition. 

To obtain further insight into the sub-gap Shiba to 
supra-gap Andreev state transition, we study the case of 
a single impurity with a classical spin S' = |S| at position 
r = 0. We use the simple gap renormalization —a^A'i5(0) 
to mimic a suppression of the gap from Aq to Aq — A' 
within a region of volume around the impurity (note 
that we require A' > 0 since there cannot be a bound 
state with energy larger than the gap if the latter is lo¬ 
cally enhanced). Inverting the Hamiltonian)^ we obtain 
an equation for the wave function at the impurity site, 
'0(0), in terms of the bound state energy E, 
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FIG. 1: The energy spectrum (black, green, and brown dots) 
and the size of the superconducting gap Ai directly under the 
impurity (red dots) for a square lattice of size 17a x 13a with 
(a) a single impurity (b) two identical impurities (Ji = J 2 ) 
as function of exchange Ji/t found numerically. Ai gets 
suppressed with increasing Ji and an effective S-S’ junction 
builds up. After the phase transition at the critical value 
Jic, Ai changes its sign giving rise to a 7r-junction. Sub-gap 
Shiba states (green dots) evolve into supra-gap Andreev states 
(brown dots) at values of Jsa where the bound state energy 
coincides with Ai. (b) For two impurities, separated by a 
distance xvija = 4, the hybridization between the two bound 
states lifts their degeneracy. There are thus two phase tran¬ 
sitions, one at Jic and one at J 2 c, where these states become 
occupied, accompanied by jumps in bound state energies and, 
hence, in Ai and A 2 . There are also two distinct Shiba to 
Andreev state transitions. The gap between the impurities 
(blue dots) is suppressed for small Ji, when the two bound 
states have finite overlap, but restores at larger Ji, see also 
Fig- El For corresponding wavefunction plots, see Fig. O The 
chosen parameters are Ao/t = 0.4 and ^,/t — 0.9. 


where a' = lypira^A', and a = Tri/pJS. We find that 
bound state energies E within the bulk gap Aq satisfying 
Eq. (fTTl) are given by 


//'o^ f (F'+ Cp'^’z + ^oTi:)(tAS' oJA'tx) ^ 

= J W- g-fp-Ag - *<“>■ 

( 10 ) 


where = p^/2m —/x, and T corresponds to the parallel 
and antiparallel orientation of the bound state spin as 
compared to the impurity spin, respectively. This equa¬ 
tion can be rewritten as 


Tx(a’ E ± oAq) -I- (o'Aq ± aE) 
\/A§ - E2 


0(0) = 0, 


( 11 ) 


E 1 — 

Aq ^ 1 -f ’ 


( 12 ) 


where r is the eigenvalue of t^,, and w = ra'ia. Plugging 
this expression for E back into Eq. (HU, we find that a 
necessary condition is |?x;| — tw = 0. Because consistency 
with BCS theory requires Aq, A' <C Ep^ and since we are 
primarily interested in the regime where |q;| is of order 
unity, we find that jal ^ |a'|. Therefore, r and ±a must 
be of the same sign, and there are two solutions with 
opposite spin like for the uniform casei^ Defining as above 
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the bound state as a Shiba state when its energy is within 
the renormalized gap, \E\ < jAp —A'|, and as an Andreev 
state when its energy is between the renormalized gap 
and the bulk gap, |Ao — A'| < \E\ < Aq, the critical 
value of the exchange coupling for the transition from a 
Shiba to an Andreev state reads 


for A' > Aq, |a| + a' > 1 (as in the numerics), and for 
A' < Aq, |a| -\-a' < 1, while otherwise the fraction under 
the square root has to be inverted. 

V. TWO IMPURITY PHYSICS 


k^'F^SA-S'l = -a' + 




A' 

2Ao - A' 




A' 

2Ao - A' 


(13) 



FIG. 2: The wavefunctions |'!/>|^ of the bound states in the case 
of a single impurity (a,c) and of two identical impurities (b,d), 
Ji = J 2 ~ J found in the tight-binding model on a lattice 
of size 17a x 17a by self-consistent numerical diagonalization. 
For two impurities, separated by a distance xvija = 4, the hy¬ 
bridization between the two bound states lifts their degener¬ 
acy. Before the phase transition with J/t = 2 (a,b), the Shiba 
bound state is well localized directly under the correspond¬ 
ing impurity. The Andreev bound state occurring after the 
phase transition at stronger values of the exchange coupling, 
J/t = 4 (c,d), is delocalized over the entire region where the 
superconducting gap is suppressed, with a maximum wave- 
function amplitude between the impurities, see (d).The cho¬ 
sen parameters are Ao/t = 0.4 and p/t = 0.9. 



FIG. 3: The cross-section of the spatial profile of the gap 
A for two impurities located on the x axis for fixed J^/t = 
2.3 and varying J\\ J\/t = 1 (yellow dotted line), Ji/t = 
2 (green dotted dashed line), Ji/t = 2.3 (red dashed line), 
J\/t = A (blue solid line). The gap under the first impurity 
is increasingly suppressed with increasing exchange coupling: 
After the phase transition Ji > Jc, tt- junctions (negative gap 
value) under the impurity arise accompanied by restoration of 
the gap between the two impurities. The gap under the second 
impurity is suppressed the most when Ji « J 2 since then the 
two bound states hybridize strongly. For wavefunction plots, 
see Fig. [5] All parameters are the same as in Fig. [T]d. 


When the impurities are further apart than the Fermi 
wavelength, ri 2 = |ri — r 2 | ^ kp^, we can expand the 
gap renormalization in orders of {kpri 2 )~^. To second 
order, we find 

c;(ri,ri,ca„)«iG“„ + G“„r«G“„ 

+ (1 + Gi2„Ti2)G2i„ (1 + , (14) 

where Tn'’ = JiSiz\l 2 x 2 - JiSizG^^^ ^ is the T- 
matrix for scattering off impurity i = 1,2. The terms 
give rise to the single impurity renormaliza¬ 
tions of the order parameter, which we denote as dA^®^ 
while the other corrections to the bare Green’s func¬ 
tion result in inter-impurity gap renormalizations 
These different contributions are depicted in Fig. |4l The 
inter-impurity renormalizations have both terms that de¬ 
pend on the relative orientation of the two spins (paral¬ 
lel for aia 2 > 0, antiparallel for aia 2 < 0), as well as 
orientation-independent contributions. All of these are of 
the order {kpri 2 )~'^. Most interesting is the orientation- 
dependent contribution evaluated at the site of one of 
the impurities. As an example, (5A2^^^(ri) reads at zero 
temperature 


(5A2^^(ri) = =F|aia2| gAo^F J (Ej 
(w2 + A2)3/2 




^F' 12 


AIA2 


{(1 - al)[{l + al)W - (1 - 


— 2(1 — Q!i)(l -I- 0 : 2 )^^ COs(2fcF?'i2)|, 


(15) 


with Ai = {1 — af)^A q + (1 + af)^w^. The upper (lower) 
sign applies for parallel (antiparallel) impurity spin align¬ 
ment. 

Besides addressing the gap renormalization at the sites 
of the impurities, it is also interesting to analyze the scal¬ 
ing of the different contributions to the gap renormaliza¬ 
tion as a function of the inter-impurity distance in their 
middle, that is at R = (ri -|- r2)/2. The single impurity 
contributions are found to scale as 


(5A«(R) ~ Ao (fcFU2)-" , (16) 


where ^ is the superconducting coherence length. This 
scaling has a simple geometrical interpretation: to lead¬ 
ing order in (kpri 2 )~^ 1, the anomalous Green’s func¬ 

tion (which determines the gap) is renormalized by the 
electron traveling from R to impurity i at r^, scattering 
there, and coming back to R. Since a trip to, and back 
from, the impurity involves a propagator proportional to 
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see Eq. dni), the gap renormalization deriving 
from the single impurity scattering processes scales as 
(fc_FCi 2 )“^. The total distance covered during these pro¬ 
cesses, precisely equal to 7 - 12 , determines the argument of 
the exponential. The leading order inter-impurity terms, 
on the other hand, are found to scale as 

5(R) ~ Ao {kpr,,)-^ , (17) 


since they stem from an electron traveling first to impu¬ 
rity one, then to impurity two, and then coming back. 
This involves three trips, and a total distance of 2ri2. 
The inter-impurity renormalization is thus sup¬ 

pressed by an additional power of kpri 2 ^ 1 in between 
the impurities as compared to its value at one of the im¬ 
purities, see Eq. (IISl). Consequently, the inter-impurity 
gap renormalization can be modeled as a func¬ 

tion with well-defined peaks close to the two impurities, 
as shown in Eig. |4l We also confirm numerically that 
the renormalizations of the superconducting gap is the 
largest directly under the impurity, while, depending on 
parameters, A(r) could even get larger than its initial 
value around the impurities, see Eig. [3l We finally note 
that the regime kpri 2 ^ 1 is qualitatively similar to the 
regime kpri 2 > 1. This is illustrated by Eig. [TJ 5 , which 
is not far from this point. 

With increasing exchange strengths, a similar transi¬ 
tion as before between Shiba and Andreev states can also 
be observed in the case of two impurities, see Fig. [TJj. 
When the two impurities are close and with comparable 
exchange strengths, the bound states overlap and their 
degeneracy gets lifted by hybridization ; which in turn 
gives rise to two separate quantum phase transitions in 
A, one at Jic and one at J 2 c-“ Similarly, we find now 
two distinct Shiba to Andreev state transitions around 
these values of exchange strengths, see Fig. [TJd. Some 
illustrative plots of the corresponding wavefunctions are 
given in Fig.[2j The spectroscopic signatures of the Shiba 
to Andreev transition are now even more drastic than 
for the single impurity case. As shown in Fig. [5 Jd, the 
bound state wavefunction of two strongly localized, but 
hybridized Shiba states is strongly suppressed between 
the impurities. In stark contrast, the wavefunction of 
the Andreev bound state spreads between the impuri¬ 
ties, and is even maximal between the impurities for the 
chosen aprameters, see Fig.[2ji. 

In the case of small impurity distances kpri 2 —>■ 0, 
where r 2 ,a;„) CIo(ri,ri,w„) = t/o(r 2 ,r 2 ,Wn). 

Keeping all orders of {kpri 2 )~^ ^ 1 hr Eq. ([7]), we find 
that the gap renormalization is in this limit given by 


i5A(ri) = 5A(r2) = (5A^^^(ri) 


Oi\—^OL—(y.\-\-OL‘2 


(18) 


where (5A*^^^(ri) is given in Eq. (|S]). Quite naturally, 
the gap renormalization resulting from two very close 
classical impurities with ai and a 2 is thus equal to the 
renormalization of a single impurity with a = ai -\- a 2 - 
Provided that J\ and J 2 have the same sign, the super¬ 
conducting gap is reduced the least if the two spins are 



FIG. 4: Sketch of the renormalization of the superconducting 
gap, 5A(r) = A(r)— Ao, in the presence of two spin-impurities 
localized at positions ri and r 2 along the x axis connecting 
the two impurities. The sketch shows the envelopes of three 
different contributions to the gap renormalization (see text). 


aligned antiparallel for short distances. As a final remark, 
we note that our non-selfconsistent analytical approach 
is also able to tackle the regime of intermediate impu¬ 
rity distances kpr 12 ^1. In this regime, we find some 
lengthy equations that simply connect the limits of small 
and large impurity distances in a smooth fashion. As an 
illustration, we refer the reader to the numerical plots in 
Fig. [Ud. These are not far from the regime kpri 2 ~ 1, 
and even treat the gap self-consistently. 


VI. CONCLUSIONS 


Analyzing two classical spins in a superconductor, we 
obtained the gap renormalization analytically for weakly 
and numerically for strongly coupled bound states. In ad¬ 
dition, we found a transition for Shiba to Andreev states 
which is accompanied by two subsequent quantum phase 
transitions as the exchange coupling is varied. 

All these predictions lead to dramatic spectral 
changes. In order to observe these effects by STM 
tecniques^id^i^SiiSl it would be preferable to control ex¬ 
perimentally the exchange interaction parametrized by a. 
This can be realized for example by considering a 2D elec¬ 
tron gas (with immersed magnetic impurities) deposited 
on top of a bulk superconductor such that the SC gap is 
now a proximity gap. We propose to bring an STM tip 
on top of a magnetic impurity in order to locally control 
the impurity density of states without affecting its spin 
and therefore to continuously monitor the parameter a. 

A number of additional directions can be envisioned 
in future work. Instead of investigating the physics as¬ 
sociated with large impurity spins, for which a classical 
treatment is appropriate, one could generalize our work 
the limit of small impurity spins, for which Kondo physics 
becomes relevant at low temperatures. With two close- 
by impurities, this could for instance lead to interesting 
multi-stage screening effects. Another promising general¬ 
ization is offered by the use of more exotic superconduc¬ 
tors, such as d-wave materials, where impurities induce 
virtual bound statesi^. 







6 


Acknowledgments 

We thank L. Glazman and S. Gangadharaiah for valu¬ 
able discussions. This work has been supported by the 


Swiss NF, NGGR QSIT, by the DFG through GRK 1621 
and SFB 1143, and by the French Agence Nationale de 
la Recherche through the ANR contracts Dymesys and 
Mistral. 


1 A. I. Rusinov, Sov. Phys. JETP 29, 1101 (1969). 

^ P. Schlottmann, Phys. Rev. B 13, 1 (1976). 

® L. Yu, Acta Phys. Sin. 21, 75 (1965). 

A. I. Rusinov, Zh. Eksp. Teor. Fiz., Pis’ma Red. 9, 146 
(1968) [JETP Lett. 9, 85 (1969)]. 

® H. Shiba, Prog. Theor. Phys. 40, 435 (1968). 

® A. Sakurai, Prog. Theor. Phys. 44, 1472 (1970). 

^ A. Yazdani, B. A. Jones, C. P. Litz, M. F. Crommie, and 
D. M. Eigler, Science 275, 1767 (1997). 

® M. E. Flatte and J. M. Byers, Phys. Rev. Lett. 78, 3761 
(1997). 

® M. 1. Salkola, A. V. Balatsky, and J. R. Schrieffer, Phys. 
Rev. B 55, 12648 (1997). 

M. E. Flatte and D. E. Reynolds, Phys. Rev. B 61, 14810 

( 2000 ). 

D. K. Morr and J. Yoon, Phys. Rev. B 73, 224511 (2006). 
A. V. Balatsky, 1. Vekhter, and J.-X. Zhu, Rev. Mod. Phys. 
78, 373 (2006). 

S.-H. Ji, T. Zhang, Y.-S. Fn, X. Chen, X.-C. Ma, J. Li, 
W.-H. Duan, J.-F. Jia, and Q.-K. Xne, Phys. Rev. Lett. 
100, 226801 (2008). 

N. Y. Yao, L. 1. Glazman, E. A. Dernier, M. D. Lnkin, and 
J. D. San, Phys. Rev. Lett. 113, 087202 (2014). 

N. Y. Yao, C. P. Moca, 1. Weymann, J. D. San, M. D. 
Lukin, E. A. Dernier, and G. Zarand, Phys. Rev. B 90, 
241108(R) (2014). 

A. A. Zyuzin and D. Loss, Phys. Rev. B 90, 125443 (2014). 
S. Hoffman, J. Klinovaja, T. Meng, and D. Loss, 
arXiv: 1503.08762 

S. Nadj-Perge, 1. K. Drozdov, B. A. Bernevig, and A. Yaz¬ 
dani, Phys. Rev. B 88, 020407(R) (2013). 


J. Klinovaja, P. Stano, A. Yazdani, and D. Loss, Phys. 
Rev. Lett. Ill, 186805 (2013). 

M. M. Vazifeh and M. Franz. Phys. Rev. Lett. Ill, 206802 
(2013). 

B. Braunecker and P. Simon, Phys. Rev. Lett. Ill, 147202 
(2013). 

S. Nakosai, Y. Tanaka, and N. Nagaosa, Phys. Rev. B 88, 
180503 (2013). 

F. Pientka, L. 1. Glazman, and F. von Oppen, Phys. Rev. 
B 88, 155420 (2013). 

K. Poyhonen, A. Weststrom, J. Rontynen, and T. Ojanen, 
Phys. Rev. B 89, 115109 (2014). 

I. Reis, D. J. J. Marchand, and M. Franz, Phys. Rev. B 
90, 085124 (2014). 

S. Nadj-Perge, 1. K. Drozdov, J. Li, H. Chen, S. Jeon, 

J. Seo, A. H. MacDonald, B. A. Bernevig, A. Yazdani, 
Science 346, 602 (2014). 

R. Pawlak, M. Kisiel, J. Klinovaja, T. Meier, S. Kawai, T. 
Glatzel, D. Loss, and E. Meyer, arXiv: 1505.06078 
A. F. Andreev, Zh. Eksp., Teor. Fiz. 49, 655 (1966) [Sov. 
Phys. JETP 22, 455 (1966)]. 

I. O. Kulik, Zh. Eksp., Teor. Fiz. 57, 1745 (1970) [Sov. 
Phys. JETP 30, 955 (1970)]. 

J. M. Rowell, Phys. Rev. Lett. 30,167 (1973). 

T. M. Klapwijk, G. E. Blonder, and M Tinkham, Physica 
B-tC 109, 1657 (1982). 

Y. Okabe and A. D. S. Nagi, Phys. Rev. B 28, 1320 (1983). 
The (non-selfconsistent) single impurity gap renormaliza¬ 
tions reported previouslyii^ are incomplete, the correct re¬ 
sult is given in Eq. ([S|l. 



